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SECTION  1 

THE  SINGER  PROGRAM 

SINGER  is  an  acronym  for  Simulation  of  Inelastic  and  Nonlinear 
Geometric  Effects  in  Reinforced  Concrete.  The  SINGER  computer  program 
addresses  the  simulation  of  a  plane  framework  assembled  of  straight 
elements,  usually  of  rectangular  cross  section.  The  user  may  model 
each  element  as  a  composite  of  longitudinal  reinforcing  bars  confined 
within  a  mass  of  concrete  by  stirrups  which,  in  turn,  have  a  protective 
cover  of  unconflned  concrete.  He  may  also  represent  wide  flange  metal 
beams.  Leaf  springs  may  be  used  to  connect  elements  eccentrically 
at  joints  or  to  model  partially  constrained  Internal  joints  and  external 
supports.  He  may  Impose  loadings  with  specified  time  histories.  He 
can  model  concentrated  and/or  distributed  forces  and  masses. 

SINGER  predicts  the  equilibrium  positions  of  the  structural  system 
either  neglecting  or  considering  inertia  effects.  It  represents  large 
changes  to  the  initial  geometry  of  the  structure  due  to  loading  and 
models  proportional,  plastic,  strain  hardening,  unloading  and  reloading 
branches  of  the  stress-strain  relations.  It  models  yielding  and  fracture 
of  materials  and  elements. 

It  characterizes  each  equilibrium  configuration  by  data  defining 
the  state  of  points  in  the  structure.  These  state  data  cite  the  displace¬ 
ment,  velocity,  acceleration,  and  fiber  strains  at  the  joints  of  the 
structure.  Additional  data  quantify  the  internal  forces  and  stresses 
at  cross-sections  where  they  may  be  excessive.  Other  output  establishes 
the  element  and  total  system  stored  and  dissipated  energy. 

SINGER  accepts  and  produces  output  in  either  English  or  System 
International  units.  (For  the  convenience  of  the  user  the  conversion 
constants  between  these  units  is  given  in  Table  1) . 
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Table  1 


CONVERSION  FACTORS  FROM  ENGLISH  SYSTEM  (EG)  OF  UNITS 
TO  INTERNATIONAL  SYSTEM  (SI)  OF  UNITS 


Unit  Multiply 

By 

To  Obtain 

Area 

in.2 

6.4516*  E-04 

2 

m 

Density  (Specific  Weight) 

lb /in. 3 

2.7144714  E+05 

N/m3 

Force 

lb 

4.4482216152605*  E+00 

N 

Length 

in. 

2.54*  E-02 

m 

Mass 

lb.s2/in. 

1.75126835  E+02 

kg 

Mass  Density 

lb.82/in.^ 

1.0686895  E+07 

kg/m3 

Mass  Moment  of  Inertia 

lb.s2in. 

1.12984829  E-01 

kg'm2 

Stress  (Pressure) 

lb /in2 

6.8947572  E+03 

N/m2 

Volume 

,  3 
in. 

1.6387064*  E-05 

m3 

*  exact,  by  National  Bureau  of  Standards  definitions. 
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SINGER  provides  special  data  manageaent  capabilities.  These  include 
logic  for  saving  data  for  a  subsequent  restart,  performing  a  restart,  and 
accumulating  a  data  retrieval  file  for  subsequent  selective  data  acquisition. 

This  guide  describes  the  Input  data  expected  and  discusses  the  output 
data  and  Its  Interpretation.  The  SINGER  Program  Document  (reference  1) 
augments  user  information  with  data  about  program  structure,  subroutine 
functions  and  flow  charts.  The  demonstration  volume  (reference  2)  describes 
demonstration  problems  and  their  results.  The  Technical  Report  (reference  3) 
details  the  mathematical  model  of  the  code. 

The  sign  convention  used  In  SINGER  for  displacements  and  forces 
depends  on  the  reference  axis  chosen.  Joint  data  (coordinates,  forces, 
displacements,  velocities,  and  accelerations)  are  referred  to  the  global 
reference  axis.  The  global  axes  must  be  chosen  so  that  the  Y-axis  is 
in  a  direction  opposite  to  gravity  (See  Figure  la).  Quantities  are 
positive  when  in  the  direction  of  the  positive  global  axes. 

Element  data  (section  properties,  distributed  loadings,  and  stress 
resultants)  are  referred  to  a  local  reference  axis.  The  local  x-axis 
coincides  with  the  reference  axis  of  the  element,  and  Its  positive 
direction  is  specified  from  the  first  listed  joint  to  the  second  listed 


1.  Barker,  R.  M.,  Melosh,  R.  J.,  and  Holzer,  S.  M. ,  "SINGER:  A 
Computer  Code  for  General  Analysis  of  Two  Dimensional  Reinforced  Concrete 
Structures",  Volume  2,  Program  Document,  Klrtland  AFB,  Albuquerque, 

New  Mexico,  Sep,  1974. 


2.  Barker,  R.  M.,  Melosh,  R.  J.,  Holzer,  S.  M.,  and  Bradshaw,  J.  C., 
"SINGER:  A  Computer  Code  for  General  Analysis  of  Two  Dimensional  Reinforced 
Concrete  Structures",  Volume  4,  Demonstration  Problems,  Klrtland  AFB, 
Albuquerque,  New  Mexico,  Sep,  1974. 

3.  Holzer,  S.  M. ,  Melosh,  R.  J.,  Barker,  R.  M. ,  and  Somers,  A.  E., 
"SINGER:  A  Computer  Code  for  General  Analysis  of  Two  Dimensional  Reinforced 
Concrete  Structures",  Volume  1,  Technical  Report,  Klrtland  AFB,  Albuquerque, 
New  Mexico,  Sep.  1974. 
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joint  (See  Figure  lb  and  Table  7) .  The  local  y  and  z-axls  coincide 
with  the  principle  axes  of  the  element.  The  local  z-axit>  ’.s  usuully 
taken  to  coincide  with  the  global  z-axls.  Quantities  are  positive 
when  in  the  direction  of  the  positive  local  axes. 

The  guide  is  organized  into  four  sections.  The  next  section 
describes  input  data  deck  form  and  card  entries.  The  third  section 
illustrates  the  output  and  discusses  its  interpretation.  The  final 
section  contains  information  on  modeling  strategy  and  suggestions  for 
controlling  problem  time  and  accuracy. 


b)  Local  (element)  coordinate  system 


Figure  1.  Coordinate  Reference  Systems 
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SECTION  2 


INPUT  DATA 

The  Input  data  consists  of  sets  of  data  for  each  problem.  Each 
i 

set  Is  composed  of  eight  data  blocks.  Each  block  consists  of  two  or 
k  more  cards  of  data.  This  section  provides  detailed  information  for 

preparing  the  input  data  blocks. 

2.1  THE  INPUT  DATA  DECK 

Figure  2  shows  how  the  problem  data  sets  form  the  Input  Data  Deck. 
Any  number  of  problems  can  be  batched  together  in  one  computer  run. 

As  long  as  all  cards  are  punched  in  the  appropriate  format,  SINGER 
will  proceed  from  problem  to  problem. 

Two  cards  terminate  the  Input  Data  Deck.  The  next  to  the  last 

« 

card  contains  the  message  END  DATA  beginning  in  Column  1  and  ending 
.  in  column  8.  This  card  activates  logic  to  indicate  that  no  further  data 

cards  are  to  be  read.  The  last  card  may  contain  zeros  or  be  left  blank. 
It  is  not  read  and  only  serves  the  purpose  of  having  all  of  the  data 
blocks  follow  the  same  pattern  of  ending  with  a  zero  card. 

Figure  3  shows  the  contents  of  each  Problem  Data  Set.  Each  deck 
includes  cards  for  each  of  the  eight  Data  Blocks.  Data  Blocks  must 
appear  in  the  order  shown  in  the  figure. 

2.2  DATA  BLOCKS 

i 

Each  data  block  contains  data  for  describing  a  particular  feature 
of  the  structural  problem  of  lnterset.  Together,  data  in  these  blocks 
either  completely  particularize  the  parameters  in  the  mathematical  model 
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Problem  Data 
Set 


Forcing  Function  Data  Block 


Function  Table  Data  Block 


Initial  Condition  Data  Block 


Point  Maas  Data  Block 


Element  Data  Block 


: 


Material  Data  Block 


Joint  Data  Block 


F  / 

'  / 


Control  Data  Block 


Figure  3.  Problem  Data  Set 
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of  the  system  or,  by  default,  evoke  values  included  in  the  computer 
code  which  represent  the  developer's  best  judgment  of  what  these  values 
should  be  to  insure  a  faithful  representation  of  a  reinforced  concrete 
framework. 

Figure  4  indicates  the  anatomy  common  to  each  Data  Block.  The  Block 
is  initiated  by  a  Block  Title  Card  which  may  contain  any  alphanumeric 
Information.  Data  on  this  card  is  printed  with  block  input  thereby 
permitting  user  Identification  of  data  in  each  block.  It  serves  no 
other  purpose.  Each  block  is  terminated  by  Zeros  Card.  This  card 
must  contain  zeros  at  least  through  column  5.  It  signals  the  end  of 
the  Block.  Numerical  Data  Cards  fall  between  these  extremes.  Even  if 
all  Numerical  Data  Cards  are  omitted,  the  Block  Title  Card  and  Zeros 
Card  must  be  in  the  Problem  Data  Set. 

Numerical  Data  Cards  are  punched  in  fixed  format  with  fields  of 
alphabetic,  integer,  or  floating  point  data  as  appropriate.  Punched 
data  may  appear  in  none  or  all  80  columns  of  the  card.  Generally,  the 
first  entries  designate  option  parameters  and  the  later  entries  provide 
problem  data. 

The  eight  tables  on  the  pages  that  follow  define  input  data  form 
and  Interpretation  for  each  of  the  eight  input  blocks  of  a  problem 
data  set.  These  tables  are  footnoted  so  they  will  be  self-contained. 
Figures  are  also  included  to  clarify  the  meaning  of  input. 

The  special  notation  below  is  used  in  these  tables  as  abbreviations: 

b  -  indicates  a  blank  (unpunched)  card  field 

d*-  indicates  the  default  interpretation-the  interpretation  made 
for  all  field  input  whose  meaning  is  not  defined  elsewhere 
in  the  input  table,  (e.g.  Table  2  indicates  that  a  "X"  punched 
in  column  16,  since  it  is  not  a  "C",  will  cause  a  STOP  if  a 
conditional  error  is  encountered  in  calculations.) 
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R  -indicates  the  card  is  required  to  be  In  the  input  deck 
(omission  of  a  R  means  the  card  may  be  omitted) 

0  -a  zero 

Card  formats  are  listed  using  FORTRAN  notation.  The  data  should  be 
right  justified  in  the  fields  specified.  The  sign  of  all  input  data 
must  be  positive  unless  otherwise  Indicated. 

Input  punched  as  zero  or  left  blank  activates  default  logic.  For 
control  information  and  option  parameters  (printout  and  simulation 
options)  default  logic  makes  a  desirable  choice  for  the  information  - 
either  what  the  program  developers  regard  as  the  most  probable  choice 
or  chat  which  will  tend  to  minimize  calculations.  In  the  case  of 
data,  default  logic  supplies,  nonzero  values,  when  necessary,  using 
standard  data.  Examples  of  default  actions  are  provided  in  the  input 


tables  which  follow. 


Card 

Cols. 

Format 

Item 

1/R 

1-80 

20A4 

Problem  Title  (to  appear  on  each  page  of 
output) 

2/R 

1-5 

15 

ISTART  -  source  of  input  data^ 

0,b  -  card  input 

l(one)  -  peripheral  unit  11 

6-10 

A2,3X 

IUNITS  -  system  of  unit  flags 

EE,d*  ■  Engllsh-ln  and  English-out 
ES  ■  English- in  and  System  Internet ional-out 
SS  -  System  Internatlonal-ln  and  System 
International-out 

SE  -  System  International-in  and  English-out 

11-15  A2,3X  ILIN  »  analysis  complexity  flag  2 

II, d*  -  inflnltesmal  joint  rotations 
(linear  analysis) 

FI  ■  finite  joint  rotations 

I  16-20  A1,4X  ISTOP  »  error  override  flag  3 

I  0,b  ■  stop  on  conditional  stop  errors'1 

C  -  continue  simulation  despite  conditional 
stop  errors 

21-25  A1,4X  IPRINT  -  print  level  flag4 

M  ■  minimum  (summaries  of  problem  description 
and  results) 

S ,d*  -  standard  (adds  input  tables  and 
solution  progress) 

D  -  detailed  (adds  unbalanced  joint  forces 
and  energies) 

E  ■  extensive  (adds  on  extended  solution 
search  printout) 

26-30  A1,4X  ISTRES  -  stress  print  flag 

N, d*  *  no  stresses  or  strains 
S  ■  stresses  and  strains 

R  ■  stress  resultants 

i  B  -  all  stresses,  strains,  and  stress 

resultants 

31-35  15  ITAPE  -  receptacle  for  continuation  data5 

O, b  -  no  continuation  data  written 
l(one)  -  peripheral  unit  11 

36-40  A1,4X  IPLOT  -  data  retrieval  file  flag 

d*  ■  no  retrieval  data  to  be  written 
R  -  write  retrieval  file  on  peripheral  unit  10 
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Table  2  (Cont . ) 


Card 

Cols. 

Format 

Item 

3/R 

1-5 

A1,4X 

I ANAL  *  Analysis  type  flag 

S,d*  ■  static  analysis 

D  -  dynamic  analysis 

6-10 

5X 

Ignored 

11-20 

E10.0 

TBEGIN  *  time  to  start  Integration,  seconds 
or  number  of  zero  load  step 

21-30 

E10.0 

THALT  *  time  to  stop  integration,  seconds  or 
number  of  last  load  step 

31-40 

E10.0 

TINK  -  approximate  time  interval  for  printing 
response  data,  seconds  of  integration 
time  (automatically  printed  every  step 
for  static  analysis) 

41-50 

E10.0 

SERR  *  maximum  tolerable  relative  error  in 
the  energy  of  unbalanced  forces® 

51-60 

E10.0 

TPR0B  -  maximum  machine  time  for  problem, 

CPU  minutes 

4/R 

1-10 

E10.0 

CA  -  ultimate  unconfined  concrete  strain 
coefficient^  (d*  *  1.0) 

11-20 

E10.0 

CB  ■  bar  buckling  end  restraint  coefficient7 
(d*  -  2.0) 

21-30 

E10.0 

CC  ■  shear  crackling  failure  curve  slope 
coefficient'  (d*  -  3.5) 

31-40 

E10.0 

CD  -  shear  cracking  axial  force  effect 
coefficient'  (d*  ■  0.0) 

41-50 

E10.0 

CE  ■  lateral  reinforcement  failure 
coefficient ^  (d*  ■  1.333) 

5/F 

1-80 

20A4 

All  zeros  or  blanks  (ends  data  block) 

Input  data  Is  either  all  from  a  card  file  or  part  from  a  card  file  and  part  from  a 
peripheral  unit.  In  the  latter  case,  the  card  file  Includes  only  the  Control  Data 
Input  Block. 


2 

Limiting  analysis  to  infinteslmal  rotations  (cos a*  1,  sina  =■  0)  reduces 
calculations  In  coordinate  transformations  and  solution  search. 

^Conditional  errors  are  Indicated  when  unusual  Input  or  calculated  data  are 
fcund.  These  may  or  may  not  be  real  errors.  For  example,  a  negative  time  to 
start  integration  Is  flagged  as  a  conditional  error  because  use  of  the 
negative  time  axis  Is  not  common. 
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4 

See  Section  3  of  this  report  for  printout  options  and  form. 

^If  continuation  data  is  written,  it  is  done  either  when  data  processing  is 
stopped  because  of  a  fatal  error  or  when  execution  continues  beyond  15  CPU 
minutes  or  if  estimate  of  calculation  time  exceeds  allowable  problem  time 
(TP ROB)  or  upon  completion  of  the  analysis. 

6  .p/2 

SERR>2  p  where  p  is  the  number  of  bits  in  the  mantissa.  This  is  the 
maximum  relative  energy  error  (deviation  from  the  minimum)  that  the  analyst 

will  accept.  The  default  value  la  about  1.  x  10  . 


^Coefficients  CA,  CB,  CC,  CD,  and  CE  modify  the  failure  criteria  limits  for 
the  expressions  developed  in  Appendix  B  of  the  Technical  Report.  Summarized 
below  in  the  English  system  of  units  are  the  expressions  which  contain  these 
coefficients.  (Section  numbers  refer  to  the  Technical  Report). 

B.3.1.1  Unconfined  Concrete  Crushing  Strain 

3+0. 002 f'  K 

efl  “  I  fVTYoob"  "  °*00351  1.0  <CA<  1.23 

where  -  unconfined  concrete  crushing  strain  and  f^  -  compressive 
cylinder  strength  of  concrete. 

B.3.1.3  Buckling  of  Longitudinal  Reinforcement 

_2E 

t  Pinned  ends  Fixed  ends 

«  ‘  16U/D)2  1.0<CB<4.0 

where  f  ■  critical  compressive  stress  in  longitudinal  bar,  E_  ■  tangent 
cr  c 

modulus  of  longitudinal  reinforcing  material,  s  ■  spacing  of  lateral 

reinforcement,  and  D  -  diameter  of  longitudinal  bar. 

B.3.2.1  Diagonal  Shear  Cracking 


v  -  1.5  /f7  +  CC •  lOOOpd  |5|  +  o  3.5  <  CC  <  5.83 

cn  c  r  1 M'  bd  -  — 


where  v  -  nominal  cracking  shear  stress  with  axial  force;  p  -  longitudinal 
cn 

steel  percentage;  d  -  effective  depth  of  element;  b  -  width  of  element; 

V,  N,  and  N  -  shear  force,  axial  force,  and  bending  moment,  respectively, 
at  critical  section;  and  a  -  coefficient  for  axial  force  effect  in  which 


tension,  ot  ■  0.025  (4  -  CD) 
compression,  acm  0.025  (4  +  2  CD) 


0  <  CD  <  1.0 


also 

lOOOpd  |jj|  <  y 


CC. 
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B.3.2.2  Yielding  of  Web  Reinforcement 

V 

where  f^  -  stress  in  web  reinforcing  bar,  ■  area  of  web  reinforcing 
bars ,  and  v  ■  nominal  shear  stress  due  to  shear  force,  V. 
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Table  3 


JOINT  DATA  BLOCK 


Card 

Cols. 

Format 

Item 

1/R 

1-80 

20A4 

Block  title  (to  lead  Joint  Data  Block  printout) 

1-5 

15 

Joint  nunber 

6-10 

5X 

Ignored 

11-20 

E10.0 

±x  coordinate  of  the  joint,  in.  or  m 

All 

21-30 

E10.0 

2 

±y  coordinate  of  the  joint,  In.  or  m 

Joint 

numerical  > 

1  31-32 

2X 

Ignored 

data 

> 

cards^- 

33 

Al 

x-displacement,  R,d*  *  restrained. 

0,b  -  unrestrained 

i 

!  34 

A1 

y-displacement,  R,d*  -  restrained. 

Restraint 
>  Codes 

i 

i 

0,b  *  unrestrained 

35 

Al 

z-rotation,  R,d*  -  restrained, 

\ 

0,b  ■  unrestrained 

4 

Last/R 

1-80 

I5,5X, 

2E10.0 

2X.3A1 

All  zeros  or  blanks  (ends  data  block) 

^Joints  i 

must  be  numbered  sequentially  beginning  from  one;  no  number  may  be 

skipped  and  the 

total  number  of  joints  equals  the  largest  joint  number. 

Joint  data  cards  need  not  be  input  In  sequential  order. 


2 

The  convention  which  must  be  followed  Is  that  the  global  Y-axls  must 
always  be  vertical  and  positive  in  the  direction  opposite  to  gravity. 
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Table  4 


MATERIAL  DATA  BLOCK1 


For 

each 

material 


Card 

Cols. 

Format 

Item 

1/R 

1-80 

20A4 

Block  title  (to  lead  Material  Data  Block  print 
out) 

'  2/R 

1-5 

A4,1X 

Material  name 

6-10 

A1,4X 

IMAT  *  material  type  flag 

U  -  unconfined  concrete 

C  ■  confined  concrete 

S,d*  -  steel 

11-20 

E10.0 

Crushing  strength  (concrete)  or  yield 
strength  (steel) *,  pai  or  N/m^,  or  f 

21-30 

E10.0 

Poisson's  Ratio,  elastic.v 

1 

• 

31-40 

E10.0 

Elastic  shear  modulus,  psi  or  N/  -,  G 

1 

1 

41-50 

E10.0 

Material  density,  lb/in^  or  N/m^,p 

i 

l 

51-60 

E10.0  3 

Cyclic  loading  curve  data 

61-70 

E10.0  J 

i 


3/R 

1-10 

A4,1X,A1, 

Card  identification^ 

4X 

N 

11-80 

7E10.0 

{  o«,  e_ ;  ...  o. 

1  1  i.  i.  M 

Coordinates  of 

.  4/R 

1-10 

A4,1X,A1, 

4 

Card  identification 

^Stress-Strain 

C 

\ 

4X 

Points 

11-80 

7E10.0 

e4*  °5»  e5*  °6»  e6’  °V 

G7 

4 

Last/R 

1-80 

A4,1X,A1, 

All  zeros  or  blanks  (ends  data  block) 

4X,  6E10.0 
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Tor  each  material  the  uniaxial  stress-strain  curve  is  modeled  by  a 
maximum  of  six  straight-line  segments  defined  by  seven  pairs  of  stress- 
strain  coordinates.  For  concrete,  no  tensile  stress  is  assumed  and 
the  compressive  values  (indicated  by  minus  signs)  are  input  as  shown 
in  Figures  5  and  6.  For  steel,  only  one  side  of  the  curve  is  given 
(Figure  7)  with  the  other  side  generated  as  a  mirror  image.  The  only 
restriction  on  the  stress-strain  curve  is  that  the  stress  must  be 
uniquely  defined  for  any  given  strain. 


2 

This  data  is  required  only  if  cards  3  and  4  are  blank  or  all  zero. 

(See  also  note  5) 

3 

This  data  particularizes  the  loading  and  unloading  curves. 

a)  For  concrete,  cols.  (51-60)  contain  the  k  of  Figures  5  and  6 
This  defines  the  slope  of  the  drop-elastic  curve,  k  must 

be  greater  than  or  equal  to  zero,  (data  in  columns  61-80 
is  Ignored). 

b)  For  steel,  cols.  (51-60,  61-70)  define  the  (a,e)  coordinates 
of  point  8  in  Figure  7.  The  unloading  curve  parallels  the 
1-2  leg.  If  reloading  occurs  with  opposite  o,  the  reloading 
curve  is  a  mirror  image  of  the  1-8-3  folded  line  of  Figure  7. 

4 

If  punched,  this  data  must  match  the  same  column  data  of  card  2. 

Thus  it  may  be  used  to  check  that  all  material  data  for  a  particular 
material  has  been  kept  together. 

^If  all  (o,e)  points  are  zero  (Indicated  by  leaving  cards  3  and  4 
blank)  the  default  stress-strain  curves  are  used.  These  are  defined 
by  the  default  data  cited  with  Figures  5,  6,  and  7.  (Table  5  ).  If 
a  steel  with  a  yield  strength  different  from  those  tabulated  is  input, 
the  stress-strain  pairs  are  automatically  Interpolated. 
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Point 

Default  Data 

Stress  a 

Strain  e 

1 

0.00 

0.00 

2 

-0.50f ' 

-0.586xl0~3 

3 

c 

-0 . 85f  * 

-1.225xl0"3 

4 

c 

-l.OOf’ 

-2.000xl0-3 

5 

c 

-0.50f' 

3+.002f ' 
c 

6 

c 

-0.35f ' 

iooo  -  r 

*  c 

7 

c 

-0.20f ' 

* 

c 


^Extrapolated  from  straight  line  through  points  4  and  5 
**Unloading  curve  slope  is  Fc  when  strain  at  which  unloading 
initiates  is  greater  than  e^.  Otherwise  the  kE  unloading 
loop  is  used. 


Figure  5.  Unconfined  Concrete  Stress-Strain  Curve 


24 


kE  **  r 

-AC  ^ 

A  4 


l  I 

1  A 

X  «.**/ 


Default  Data 


Point 

Stress .o 

Strain  ,e 

1 

0.00 

0.00 

2 

-0.50f ' 

C 

-0.586  x  10"3 

3 

-0.85f ' 

-1.225  x  10“3 

c 

4 

-i-oofl 

c 

-2.000  x  10‘3 

5 

-<f;  +  Afc> 

-(.002  +  Ae  ) 
c 

5.5 

-0.50f ' 

-(3+.002fg)  3 

c 

(f’-lOOO.)  4* 

6 

-0.20f * 

C 

* 

7 

-0.20f  ’ 

C 

e?  -  0.3 

ii  'r. 

S 


Af  -  xP"f  " 
c  4  r  y 

Aec  -  O.HpM/Ts 

p"  -  ratio  of  volume  of  lateral  reinforcement  to  confined  concrete 
1^"  -  yield  stress  of  lateral  reinforcement 
b"  »  width  of  confined  concrete 
s  -  spacing  of  lateral  reinforcement 

^Extrapolated  from  a  straight  line  through  points  5  and  5.5 

**Unloading  curve  slope  Is  E£  when  strain  at  which  unloading  initiates  is 
greater  than  e^.  Otherwise  the  kE£  unloading  loop  is  used. 

Figure  6.  Confined  Concrete  Stress-Strain  Curve 


a 


Figure  7.  Steel  Stress-Strain  Curve 
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♦♦Unloading  curve  point  (point  8  of  Figure  7 ) .  Specific  values  are  obtained  from  this  table  on  default  using  linear 
Interpolation.  If  extrapolation  is  required,  an  error  message  is  produced. 


Table  6 


ELEMENT  DATA  BLOCK1 


Card  Cols .  Format 

1  1-20  20A4 

l/R3 

2 

3,4, . . . L(max  L-12) 

(L-no.  of  longitudinal  bar  groups) 

L+l,L+2,...M(M<L+6) 

(M  -  no.  of  lateral  bar  groups) 

M+l 

M+2 

Last/R  1-5  15 


Item 

Block  title  (to  lead  Element  Data  Block 
printout) 

Element  parameter  card  1  (see  Table 
7  > 

Concrete  data  card  (see  Table  8  ) 

Longitudinal  reinforcement  card 
(see  Table  9  ) 

Lateral  reinforcement  card^  (see 
Table  10  ) 

Hide  flange  reinforcement  card  (see 
Table  11  ) 

Leaf  spring  flexibility  card  (see 
Table  12) 


All  zeros  or  blanks  (ends  data  block) 


Groups  of  cards  for  each  element  can  be  input  in  any  order.  All  cards  in 
the  group  have  the  same  format. 


2 

Except  as  noted  in  footnote  4,  each  element  parameter  card  may  be  followed 
by  any  of  the  other  card  types  listed.  At  least  one  such  card  must  be 
Included  and  cards,  if  given,  must  follow  the  sequence  given. 


3 

This  card  is  omitted  for  a  vide  flange  or  leaf  spring  element. 

4 

Lateral  reinforcement  cards  cannot  be  introduced  unless  preceded  by  one  or 
more  longitudinal  reinforcement  cards. 
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Table  7 


ELEMENT  PARAMETER  CARD 


Cols. 

Format 

Item 

1-5 

A4,1X 

Blank 

6-10 

15 

Element  number^ 

11-15 

15 

Joint  number  for  one  end  of  the  member,  or 
the  "F"  (first)  end. 

16-20 

15 

Joint  number  for  the  other  end  of  the  member, 
or  the  "S"  (second)  end 

21 

A1 

1ACT  -  element  action  flag 

L,  -  restricted  linear  material  behavior  (Hooke's 
law) 

C,d*  -  changeable,  initially  Hookean  but  may 
become  nonlinear  due  to  loading 

N  -  non-linear  material  behavior 

22 

A1 

ISHEAR  -  shear  model  flag 

N,d*  ■  don't  check 

S  -  check  using  linear  behavior  range 

23 

A1 

2 

I BOND  ■  steel  bond  check  flag  ACI  318-71 

N,d*  »  don't  check 

H  -  check  using  ACI  318-71  specs 

24-30 

Al,  11,  A4,  IX 

Read  and  ignored 

31-40 

E10.0 

±h  *  distance  to  reference  axis  from  top  of 

C  3 

element  ,  in.  or  m. 

41-50 

E10.0 

3 

±d  -  distance  to  centroid  of  bottom  rebar  group 

2 

from  top  of  element  ,  In.  or  m. 

51-60 

E10.0 

3 

±d'  ■  distance  to  centroid  of  top  rebar  group 

2 

from  top  of  element  ,  In.  or  m. 

61-70 

E10.0 

Xp  -  distance  from"F"end  to  section  used  in  failure 
check**,  in.  or  m. 

71-80 

E10.0 

Xc  ■  distance  from  "S"  end  to  section  used  In 

3 

failure  check  end  ,  in.  or  m. 

*If  no  element  number  la  specified  by  the  user,  elements  are  assigned  numbers 
corresponding  to  their  order  In  Input. 


29 


This  check  is  a  comparison  of  the  effective  length  of  the  reinforcement 
group  with  the  development  length  required  to  reach  a  bar  stress  of 
f  (Sec.  12.5,  ACI  318-71). 

3 

Figure  8,  below,  clarifies  the  meaning  of  these  input  data.  The  X 
distances  are  from  the  joint  to  the  section  (usually  at  face  of  supports) 
where  flexure  and  shear  failure  checks  are  made.  The  Xt  values  do  not 

effect  the  prediction  of  Joint  displacements.  The  reference  axis  position 
is  arbitrary  and  immaterial  to  response  prediction.  Either  end  may  be  the 
first  or  second  end. 


Figure  8.  Element  Gross  Geometry 
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Table  8 


CONCRETE  DATA  CARD 


Cols. 

Format 

Item 

1-5 

A4.1X 

Unconfined  concrete  material  name^ 

6-10 

15 

2 

Element  number 

11-20 

215 

Read  and  ignored 

21-25 

4A1.I1 

Read  and  ignored 

26-30 

A4,1X 

Confined  concrete  material  name^ 

31-40 

E10.0 

3 

h  -  element  height  ,  in.  or  m. 

41-50 

E10.0 

3 

b  -  element  width  ,  in.  or  m. 

51-60 

E10.0 

d"  -  confined  depth  of  element  ,  in.  or  m. 

61-70 

E10.0 

3 

b"  -  confined  width  of  element  ,  in.  or  m. 

71-80 

E10.0 

Read  and  ignored 

Rincon fined  concrete 

is  that  not  bound  by  stirrups.  The  confined  concrete 

goes  within  the  etlrrups.  Material  data  must  be  included  in  the  Material 
Data  Block  for  each  concrete  name. 

2 

If  no  element  number  ia  specified  by  the  user,  this  number  will  be  that  of 
the  preceding  parameter  card  (see  Table  7  ) . 

3 

Figure  9,  on  the  next  page,  clarifies  the  meaning  of  these  input  data. 
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Note  relationship  between  h  and  the  element's  y-axis.  The  element's  x-axis 
points  into  paper. 

When  setting  up  a  composite  member  with  an  external  cover  plate,  the 
element  height  (h)  is  to  be  given  as  the  height  of  the  concrete  section  (i.e, 
it  does  not  include  the  tnickness  of  the  plate).  The  cover  plate  area  and 
distance  from  the  reference  axis  is  specified  by  way  of  a  longitudinal  rein¬ 
forcement  card  (Table  9  ) . 


Figure  9,  Element  Cross  Section  Dimensions 
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Table  9 


LONGITUDINAL  REINFORCEMENT  CARD1 


Cols. 

Format 

Item 

1-5 

A4,1X 

"BARS"  -  Alphanumeric  to  indicate  that  longitud¬ 
inal  reinforcement  is  to  be  input 

6-10 

15 

2 

Reinforcement  group  number 

11-15 

15 

4 

Number  of  bars 

16-20 

15 

Bar  size  number1’^ 

21 

A1 

Reinforcement  continuity  code  at  the  "F"  end  of 
the  current  reinforcement  group;  0-reinforcement 
terminated;  1-rebar  bent  or  continued 

22 

A1 

Reinforcement  continuity  code  at  the  "S"  end  of 
the  current  reinforcement  group,  0  or  1 

23-25 

2A1.I1 

Read  and  ignored 

26-30 

A4.1X 

Reinforcement  material  name 

31-40 

E10.0 

4  2  2 

Area  of  reinforcement  group  ,  in.  or  m. 

41-50 

E10.0 

y^  -  distance  to  steel  from  reference  axis5, 
in.  or  m. 

51-60 

E10.0 

X  -  distance  from  ’V"  end  of  member  to  beginning 
of  reinforcement  group^,  in.  or  m. 

61-70 

E10.0 

X  -  distance  from  end  of  reinforcement  to  "S"  end 
s  element,  in.  or  m. 

71-80 

E10.0 

Read  and  ignored 

up  to  ten  longitudinal  reinforcement  groups  are  allowed  on  each  element. 

One  Longitudinal  Reinforcement  Card  is  read  for  each  group. 

2 

The  user  may  number  the  reinforcement  groups.  The  numbers  must  be  specified 
without  skipping  numbers  and  starting  from  1.  If  the  reinforcement  group 
number  is  not  specified,  groups  are  assigned  numbers  corresponding  to 
their  order  in  input. 


3 

If  metric  bar  sizes  are  being  used,  the  bar  size  number  should  be  the 
diameter  of  the  bar  in  millimeters. 

4 

Two  options  are  provided  for  calculating  the  reinforcement  group  area. 

The  user  may  input  the  area  of  the  group  directly  or  the  number  of  bars 
and  the  bar  size  number.  If  the  user  inputs  the  number  of  bars  and 
the  bar  size  number,  subroutine  TABL  will  calculate  the  area  and  the 
bar  size  diameter,  if  the  bar  size  is  specified,  it  will  (optionally) 
determine  if  bond  is  sufficient,  and  if  reinforcement  buckling  is  likely 
to  occur.  If  only  the  group  area  is  input,  these  calculations  are  skipped. 


Figure  10.  Longitudinal  Reinforcement  Dimensions 
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Table  10 


LATERAL  REINFORCEMENT  CARD1 


Cols. 

Format 

Item 

1-5 

A4,1X 

"STIR"  ■  stirrup  group  or  "TIES"  ■  tie 
group,  identification2 

6-10 

15 

3 

Lateral  reinforcement  group  number 

11-15 

15 

4 

Bar  size  nuaber  of  lateral  reinforcement 

16-20 

15 

Number  of  spaces  in  each  stirrup  or  tie  group  (See  Fig.  11) 

21-24 

4A1 

Read  and  ignored. 

25 

11 

Number  of  legs  for  individual  stirrup  or  tie 

26-30 

A4.1X 

Lateral  reinforcement  material  name 

31-40 

E10.0 

Ay  «  total,  vertical  area  of  lateral  reinforcement 
group3,  in. 2  or  m2. 

41-50 

E10.0 

S  «  spacing  of  stirrups  in  this  group,  or  spacing 
of  ties,  in. ,  or  m. 

51-60 

E10.0 

X  ■  distance  fron  "F"  end  of  element  to  beginning 
the  current  tie  or  stirrups  spacing,  in.  or  m. 

61-70 

E10.0 

5  3  3 

V  ■  volume  of  steel  per  stirrup  ,  in.  or  m.  . 

0 

71-80 

E10.0 

Read  and  ignored 

^rom  1  to  6  stirrup  or  tie  groups  may  be  specified  for  lateral  reinforcement.  One 
lateral  reinforcement  card  is  read  for  each  such  group. 

2 

As  far  as  the  analysis  Is  concerned,  there  Is  no  distinction  between  lateral  rein¬ 
forcement  designated  as  "STIR"  or  "TIES”.  Two  terms  are  provided  to  permit  the 
user  to  use  conventional  terms  -  stirrups  for  beams  and  ties  for  columns. 

3 

The  stirrup  or  tie  groups  may  be  specified  by  the  user.  If  not,  groups  are  assigned 
numbers  corresponding  to  their  order  in  input. 

4 

If  metric  bar  sizes  are  being  used,  the  bar  size  number  should  be  the  bar  diameter 
in  millimeters. 


^ If  the  user  wishes  he  may  input  the  area,  A  ,  directly  and  not  input 
the  bar  size  and  number  of  legs.  If  both  entries  are  used,  the  program 
determines  the  area  and  volume  of  stirrups.  If  only  the  area  is  given, 
the  user  must  also  supply,  Vs,  the  volume  of  steel  per  stirrup. 


Figure  11.  Lateral  Reinforcement  Dimensions 
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Table  11 


WIDE  FLANGE  REINFORCEMENT  CARD 


Cols. 

Format 

Item 

1-5 

A4,1X 

"WFST"  "  wide  flange  card  Identification 

6-10 

15 

Element  number^ 

11-15 

15 

2 

Joint  number  of  the  "F"  end  of  the  element 

16-20 

15 

2 

Joint  number  of  the  "S"  end  of  the  element 

21 

A1 

? 

IACT  -  element  action  flag~;  L,  C(d*>  ,  or  N  (see  page  29) 

22-25 

3A1.I1 

Read  and  ignored 

26-30 

A4.1X 

Wide  flange  material  name 

31-40 

E10.0 

h  ■  Distance  from  top  surface  to  reference  axis, 
in.  or  m. 

41-50 

E10.0 

h  »  Depth  of  element,  in.  or  m. 

51-60 

E10.0 

t  ■  Thickness  of  element  web,  in.  or  m. 
w 

61-70 

E10.0 

b  -  Width  of  element,  in.  or  m. 

71-80 

E10.0 

t.  -  Thickness  of  element  flange,  in.  or  m. 

If  no  element  number  la  specified  by  the  user,  this  number  will  be  assigned 
as  that  of  the  preceding  element  parameter  card  (see  Table  7  )  or  If  this  Is 
zero,  the  next  number  In  sequence. 


2 

If  joint  numbers  and  IACT  are  omitted,  the  values  will  be  taken  from  the  previous 
element  parameter  card  and  the  wide  flange  beam  will  be  regarded  as  the  steel 
reinforcement  of  a  reinforced  concrete  element. 


h 


Figure  12. 


37 


Table  12 


LEAF  SPRING  FLEXIBILITY  CARD1 


Col8. 

Format 

Item 

1-5 

A4,1X 

"LEAF1 

"  ■  leaf  spring  card  identification 

6-10 

15 

2 

Leaf  spring  number 

11-15 

15 

3 

Joint  number  of  F  end  of  spring 

16-20 

15 

3 

Joint  number  of  S  end  of  spring 

21-30 

4A1,  11,  A4,  IX 

Read  and  ignored 

31-40 

E10.0 

*11  “ 

lateral  deflection  for  unit  tip  lateral 
load^,  in. /lb.  or  m./N 

41-50 

E10.0 

A12 

z  axis  rotation  for  unit  tip  lateral 
load^,  rad/lb.  or  rad/N 

51-60 

E10.0 

*22  " 

4 

z  axis  rotation  for  unit  tip  z  moment  , 
rad. /in. lb.  or  rad./m.N 

61-70 

E10.0 

*33“ 

element  x  axis  elongation  for  unit  tip 
axial  load,  in. /lb.  or  m./N 

71-80 

E10.0 

Leaf  spring  storable  energy  limit ^ »  in-lb  or  ' 

1Deflnes  load-deflection  coefficients  for  a  linear  leaf  spring.  Using  a  leaf 
spring  to  connect  two  elements  provides  for  Implying  any  degree  of  Interelement 
continuity. 

2 

Leaf  springs  have  a  numbering  system  which  Is  separate  from  those  assigned  to 
other  elements.  This  numbering  begins  with  one  and  then  proceeds  In  sequence. 


3 

Since  two  joint  numbers  must  appear  on  this  card,  it  can  be  put  anywhere  In 
the  element  data  block. 

4 

Flexibility  data  is  interpretated  as  influence  coefficients  for  forces  and 
moments  as  illustrated  in  Figure  13,  on  the  next  page.  Note  that  the  flexibility 
data  Is  defined  relative  to  the  element  local  coordinates  (as  defined  in  Table  16, 
footnote  2).  (At  most,  this  axis  direction  only  effects  the  sign  of  A^*) 

^When  leaf  spring  behavior  involves  strain  energy  equal  or  greater  than  this  limit, 
a  warning  wssage  is  printed.  An  extremely  large  limit  inhibits  the  warning. 
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Unit  Force 

Fig.  13b  Fig.  13d 


Figure  13 ,  Leaf  Spring  Flexibility  Data 


Table  13 


POINT  MASS  DATA  BLOCK1 


Card 

Cols. 

Format 

Item 

1  /R 

1-80 

20A4 

Block  title  (to  lead  Point  Mass  Data  Block) 

f 

1-5 

15 

Joint  number 

One  card  for 
each 

6-10 

5X 

Not  read 

concentrated  v 

mass  (minimum 
none) 

11-40 

i 

1 

3E10.0 

M  -  mass  parameter  for  motion  in  the 
x  x-direction  ,  lb.s?/in. ,or  kg. 

M  -  mass  parameter  for  motion  in  the 
y  y-direction  ,  lb.s.2/in.,  or  kg. 

I  *  mass  moment  of  inertia  for  rotation 

2  about  z-axis,lb.s2in.  or  kg.m2. 

Last  /R 

1-5 

15,  5X 
3E10.0 

All  zeros  or  blanks  (ends  data  block) 

The  program  generates  lumped  masses  due  to  element  weight  when  the  material 
density  is  non-zero,  except  for  leaf  spring  elements.  If  their  mass  is  to  be 
represented,  it  is  Introduced  as  point  mass  data. 


40 


Table  14 


INITIAL  CONDITION  DATA  BLOCK1*2 


Card 

Cols. 

Format 

Item 

1/R 

1-80 

A 

Block  title  (to  lead  Initial  Condition  Data  Block) 

✓ 

1-5 

15 

Joint  number 

Cards  2 

through 

s 

next  to 

6-10 

A1,4X 

type  of  initial  condition 

D,d*  -  displacement3,  in.  or  m. 

V  ■  velocity, in/sec  or  m./s 

A  ■  acceleration  in. /sec*  or  m/sec 

J  ■  jerk  in/sec3  or  m/sec3 

F  -  force*»5  lb  or  N. 

last 

11-20 

E10.0 

x -direction  component 

21-30 

E10.0 

y-direction  component 

31-40 

E10.0 

z-axls  rotation  component 

Last/R 

1-5 

15, Al, 
4X, 

3E10.0 

All  zeros  or  blanks  (ends  data  block) 

'If  no  Initial  conditions  are  specified,  (Cards  2  through  next  to  lest  are  omitted), 
a  static  analysis  Is  performed. 


2 

Initial  conditions  for  all  degrees  of  freedom  at  joint  are  either  "given"  or 
"unspecified".  They  are  "given"  by  Including  an  initial  condition  card  for 
the  joint.  Conditions  for  all  degrees  of  freedom  at  the  joint  are  taken  from 
the  card.  They  are  implied  to  be  "unspecified"  if  no  Initial  condition  card 
references  the  joint. 


3 

A  one  g  gravitational  acceleration  Is  assumed  to  act  in  the  negative  y  direction. 
This  Induces  inertia  forces  for  all  structural  material  with  a  non-zero  density  In 
the  material  table. 

4 

When  Initial  displacements  or  forces  are  given,  SINGER  does  a  time  zero  static 
analysis.  vor  initial  displacements,  this  produces  associated  initial  stresses. 
For  Initial  forces,  this  yields  an  initial  state  of  displacement.  When  the 
integration  starts,  these  forces  are  removed.  (The  option  of  giving  initial 
forces  relieves  the  user  of  calculating  an  initial  displacement  state). 

^If  more  than  one  set  of  forces  are  given  for  a  particular  joint,  forces  are 
accumulated  to  develop  a  set  of  total  joint  forces. 
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Table  15 


FUNCTION  TABLE  DATA  BLOCK1 


Card 

Cols. 

Format 

Item 

1/R 

1-80 

20A4 

Block  title  (to  lead  Function  Table  Data  print 
out) 

2/R 

1-5 

14, IX 

2 

Function  number 

6-10 

5X 

Ignored 

11-70 

6E10.0 

^1*  ^1*  fc2»  *2*  * ** 

3 

function  coordinates 

3/R 

1-5 

15 

2 

Function  number 

6-10 

5X 

Ignored 

11-70 

6E10.0 

•  •  • »  b  ,  F 
n  n 

3 

function  coordinates 

Last  /R 

1-5 

14, 6X 
6E10.0 

All  zeros  or  blanks 

(ends  data  block) 

Hhe  function  tables  define  reference  functions  for  describing  tine  or 
load  step  variations  of  joint  forcing  functions. 


2 

The  functions  must  be  numbered  sequentially  beginning  at  one;  no  number  may 
be  skipped  or  omitted;  the  total  number  of  functions  equals  the  largest 
function  number.  The  function  data  cards  must  be  ordered  sequentially 
by  function  number 
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3 

Each  function  is  defined  by  straight-line  segments  connecting  prescribed 
points,  as  shown  in  Figure  14.  There  is  no  explicit  limit  on  the 
maiber  of  points. 

(There  is  an  illicit  limit  on  the  number  of  points.  The  function  tables 
data  is  stored  dynamically  in  data  arrays  along  with  other  program  data. 

When  the  array  size  is  exceeded  an  error  message  is  printed.) 

The  function  must  be  single-valued  in  time  or  load  step;  "backwards"  slopes 
are  not  permitted.  For  vertical  slopes,  the  average  value  of  the  two  end 
points  is  used. 


Figure  14.  Function  Form 
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1 

I 


Table  16 

FORCING  FUNCTION  DATA  BLOCK1 


Card 

Cols. 

Format 

Item 

1/R 

1-80 

20A4 

Block  title  (to  lead  Forcing  Function  Data  Block) 

f 

1-5 

15 

2 

Joint  nunber  loaded  joint  ,p 

6-10 

15 

Joint  number  of  reference  joint  •  Q 

11-15 

A1,4X 

P,d*  -  distributed  loading  intensity  at  point  p 
(normal  to  element  axis),  lb /in^  or  N  /M 

X  ■  x-direction  concentrated  force*,  lbs  or  N 

Y  m  y-direction  concentrated  force  ,  lbs  or  N 

Z  ■  z-direction  concentrated  moment,  lbs  or  N 

each  force 

16-19 

14 

Referenced  function  nunber 

component 

■ 

20 

A1 

K,d*  -  constant  function  of  magnitude  LA 

S  ■  sinusoidal  function 

C  ■  cosinusoidal  function 

21-60 

i 

j 

s 

4E10.0 

SF,  LA,  TA,  TP3 

SF  ■  scale  factor  (d*  »  1.0) 

LA  ■  load  addition  constant*  or 

TA  ■  time  addition  constant,  s. 

TP  ■  time  period  for  sinusoidal  or  cosinusoidal 
functions  (not  used  otherwise),  s. 

Last  /R 

1-5 

215,  A1.4X, 
I4,A1,4E10.0 

All  zeros  or  blanks  (ends  data  block) 

1The  joint  forcing  function  data  defines  time  or  load-step  varying  forcing  functions 
applied  at  the  joints  or  on  elements. 


When  a  concentrated  load  is  desired,  only  the  loaded  joint  number  is  needed,  so  the 
reference  joint  number  is  Ignored.  Specification  of  the  reference  joint  number 
Indicates  a  linearlly  distributed  load  is  Involved  with  the  magnitude  of  the 
pressure,  per  unit  length,  specified  for  the  loaded  joint  Pand  varying  linearlly  to 
a  zero  value  at  the  reference  joint  Q.  This  pressure  is  assumed  to  act  normal  to 
the  element  reference  axis.  A  positive  pressure  is  outward  to  the  left  when 
proceeding  from  P  to  Q.  The  relation  of  this  direction  to  the  P-Q  element 
ends  and  the  global  axis  is  illustrated  in  Figure  15 . 
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1 


p 


Figure  15  .  Positive  Direction  of  Distributed  Pressure 


Function  Scaling  Parameters 


3 

The  time  variation  is  described  by  a  reference  function  as  suggested  by 
Figure  16.  The  reference  function  obtained  from  the  function  tables,  may  be 
scaled  by  a  multiplication  factor  (SF)  .shifted  in  magnitude  by  an  addition 
factor  (LA),  or  shifted  in  time  by  a  time  factor  (TA).  Multiple  forcing 
functions  at  a  joint  are  superimposed. 
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2.3  IMPLICIT  MODELING  FEATURES 


SINGER  contains  special  logic  corresponding  to  three  of  the  input 
blocks.  This  logic  makes  implicit  in  the  input  specifications  the  ability 
for  refining  the  finite  element  model,  introducing  joint  eccentricities  at 
connections,  and/or  representing  linearly  varying  pressure  loadings  only 
by  adding  to  the  input  data.  The  next  paragraphs  clarify  how  the  user 
can  prepare  input  to  exploit  these  features  of  the  input  specification. 

2.3.1  Addition  of  Interior  Joints 

The  input  specification  provides  for  a  complete  description  of  the 
structural  configuration.  Occasionally,  however,  it  is  necessary  to  rep¬ 
resent  the  system  with  a  more  refined  mathematical  model.  This  may  be 
prompted  by  SINGER  output  or  the  user's  desire  for  greater  accuracy,  for 
example.  This  is  accomplished  by  adding  interior  joints  to  the  joint  data 
block. 

An  interior  joint  is  a  joint  that  is  on  the  axes  of  an  element 
between  the  F  and  S  terminating  joints.  In  input  data,  it  is  lndls- 
tlnqulshable  from  other  joints  in  the  input  data  block.  It  is  detected 
by  discovering  that  the  joint  lies  on  a  line  between  a  pair  of  element 
joints. 

When  an  interior  joint  is  detected,  it  causes  replacement  of  the 
pertinent  element  data  by  data  for  two  elements  with  the  same  aggregate 
span  and  characteristics.  Once  this  has  been  done,  the  interior  joint 
loses  its  identity.  Because  of  this,  the  interior  joint  logic  has  no 
limit  on  the  number  of  interior  joints  the  user  may  require  between  any 
two  joints  of  the  configuration  descriptive  data. 
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Loss  of  the  identity  of  interior  joints  also  means  that  output  does 
not  make  this  distinction.  Joint  displacements  for  interior  joints  appear 
with  the  interior  joint  number.  Element  behavior  is  specified  for  both 
parts  of  a  subdivided  element  just  as  if  the  user  had  prepared  the  input 
in  the  form  of  two  separate  sets  of  element  data. 

Note  that  an  interior  joint  is  only  effective  in  generating  additional 
element  data  and  associated  printouts.  It  does  not  change  the  mass  or 
joint  loading  data,  even  for  distributed  loads.  If  it  is  necessary  to 
refine  the  materials,  mass  or  loading,  the  user  must  revise  the  input 
data. 


2.3.2  Joint  Eccentricities 

The  ability  to  represent  joint  eccentricities  is  implicit  in  the 
input  specifications  for  element  data.  Eccentricities  are  defined  using 
leaf  spring  input. 

For  example,  suppose  it  is  desirable  to  model  an  eccentric  connection 
between  two  elements:  one  with  F  and  S  joints  1,  2;  the  other  with  joint 
3,  4.  Then,  if  the  eccentric  connection  is  between  joints  2  and  3,  a  leaf 
spring  with  zero  flexibility  coefficients  is  added  to  the  input.  If  the 
eccentricity  does  not  exist,  of  course,  joint  2  and  3  would  be  assigned 
the  same  joint  number  by  the  analyst. 

Note  that,  unlike  linear  analysis  simulators,  the  representation  of 
eccentricities  is  not  "exact".  Accordingly,  output  describing  response 
will  imply  as  ouch  flexibility  in  the  leaf  spring  as  the  user  permits  by 
his  specification  of  tolerance  on  energy  error  in  the  Control  Data  Block 
Input  Data. 
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2.3.3  Linearly  Varying  Pressure  Loadings 

A  uniformly  distributed  load  Is  represented  as  the  superposition  of 
two  distributed  loads  in  the  Forcing  Function  Data  Block.  The  two  loads 
must  involve  the  same  joint  pairs  for  the  loaded  and  reference  joint. 

Suppose,  for  example,  element  1-2  is  loaded  with  a  uniform  pressure 
of  intensity  q.  One  joint  loading  card  vrould  indicate  joint  1  as  the 
loaded  joint  and  2  as  the  reference.  The  lead  Intensity  would  be  q. 

The  other  joint  loading  card  would  Involve  joint  2  as  the  loaded  joint  an 
1  as  the  reference.  In  this  case,  because  the  loading  is  outward  to  the 
right  going  from  the  loaded  joint  to  the  reference  joint,  the  load 
Intensity  would  be  -q. 

As  illustrated  in  Figure  17,  the  superposition  of  these  two  loadings 
represents  the  uniformly  distributed  load. 

The  superposition  permits  representing  any  llnearllly  varying  load 
between  two  joints.  A  given  linearly  distributed  load  can  be  fractionated 
into  as  many  superlmposable  components  as  desirable. 

Note  that  the  distributed  loading  is  not  represented  exactly  in  the 
mathematical  model,  just  as  the  distributed  mass  of  the  structure  is 
approximated.  Exact  representation  would  require  evaluating  equilvalent 
joints  loads  for  each  deformation  state.  Accordingly,  SINGER  logic  may 
caution  the  user  to  Introduce  interior  joints  or  Improve  the  joint  loading 
input  resolution  for  a  more  accurate  simulation. 
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Joint  1  First  Loading  Joint  2 

(P)  (Q) 


^rrrrm™- 

Joint  1  Second  Loading  Joint  2 

(Q)  (P) 


Joint  1  Total  Loading  Joint  2 


Figure  17.  Superposition  of  Distributed  Loadings 
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SECTION  3 


OUTPUT  DATA 

Output  data  Includes  both  printed  data  and  data  stored  on  peripheral 
storage  units.  Printed  results  define  the  tine  history  of  the  structural 
configuration  changes  and  data  transferred  to  retrieval  files  on  peripheral 
storage  units  facilitates  plotting  structural  behavior  and  restarting  the 
analysis.  This  section  describes  and  Illustrates  the  form  and  interpreta¬ 
tion  of  this  output. 

3.1  PRINTED  OUTPUT 

Printer  output  consists  of  a  definition  of  the  problem  statement, 
available  after  all  Input  haa  been  read  and  checked;  data  describing  the 
state  of  the  system,  printed  upon  completion  of  calculations  for  each  time 
period;  and  diagnostic  messages,  printed  when  anomalies  arise  during  the 
simulation. 

The  amount  of  printed  output  depends  on  the  options  elected  and 
specified  In  the  Control  Data  Block.  The  amount  depends  on  the  print¬ 
out  increment  level  and  the  dimensional  units  Involved  in  input  and  output. 

Table  17  lists  the  Increments  of  printed  output  for  each  level 
selected  and  for  each  group  of  output.  For  a  given  level  of  output,  all 
data  at  that  and  all  lower  levels  Is  printed.  Thus,  as  the  table  shows, 
at  the  detailed  option  level,  Input  data  printouts  include  the  table 
of  the  input  deta  for  each  input  data  block  as  well  as  the  summary 
description  of  the  problem  under  study. 

When  the  dimensional  units  for  Input  and  output  are  mixed  (i.e.  one  is 
In  English  units  and  the  other  Is  System  International  units)  the  amount  of 
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Table  17 


PRINTOUT  INCREMENTS  FOR  EACH  OPTION  LEVEL 


Option 

Level 


1.  Minimum 

2 .  Standard 

3.  Detailed 

4 .  Extensive 


Input  Data 
Increment 

Problem  descrip¬ 
tion  8  usury 


Input  tables 
of  Block  Data 


State  Data* 
Increment 

Joint  response 
Element  energies 
Accuracy  susmary 

Solution  progress 

summary 

Unbalanced  joint 
forces  and  energies 

Detailed  solution 
search  progress 


Message 

Increment 

Fatal  and  critical 

mesaages 


Potentially  critical 

messages 


Element  Internal  forces  and/or  stresses  and  strains  are  printed  If  required 
by  the  ISTRES  flag  of  Block  1  Input  data,  regardless  of  the  printout  option 
level. 
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printed  output  is  increased.  The  mix  incites  printing  of  the  input 
data  in  both  sets  of  units  to  provide  a  copy  of  input  data  in  the 
same  dimensional  units  as  the  output. 

Each  page  of  printed  output  has  a  heading,  subtitling  and  problem 
data  or  comments.  The  page  heading  consists  of  a  centered  copy  of  the 
first  card  of  input  from  Control  Data  Block  and  the  page  number.  Titling 
designates  the  group  of  output.  The  problem  data  includes  tabulated 
numbers,  and  names  with  column  labels  defining  dimensional  units,  and 
relevant  diagnostic  messages. 

There  are  a  maximum  of  120  characters  per  line  of  printout.  Type¬ 
written  notes  do  not  appear  on  the  computer  printout  but  have  been  added 
to  provide  a  self-contained  interpretation  of  page  data. 

Actual  SINGER  printouts  may  differ  from  these  tables  because  the  ones 
given  here  were  produced  by  the  original  version  of  the  code  on  the  IBM  370. 
This  version  prints  a  maximum  of  66  lines  per  page,  uses  the  EBCDIC 
character  set,  and  prints  from  double-precision  floating  point  data. 

3.1.1  Printed  Input  Data 

Printed  input  consists  of  seven  groups  of  pages;  one  group  corresponding 
to  each  of  the  last  seven  input  blocks.  An  eighth  group  of  pages  defines 
the  finite  element  analysis  articulation.  The  last  group  of  pages  of 
input  printouts  provides  a  summarizing  description  of  the  structural 
configuration  under  study. 

The  example  problem  of  Figure  18  is  used  to  describe  the  nature  of 
the  printed  input  data  developed  by  SINGER.  Tables  18  through  24  illustrate 
the  listings  which  correspond  to  the  last  seven  input  blocks.  These  data 
completely  describe  the  structural  configuration  -  geometry,  materials, 
and  boundary  conditions. 
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Table  18.  Joint  Coordinate  Data 
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The  subtitling  for  each  group  contains  the  block  title  prepared 
by  the  user  and  Included  as  the  first  card  in  the  input  block.  This 
title  is  reproduced  just  as  it  appears  on  the  input  card.  It  is  printed, 
left  justified.  (If  the  user  wants  this  title  centered,  he  oust  center 
it  about  column  60  on  the  card  input) . 

Subtitling  also  Includes  the  labels  which  lead  each  column  of 
tabular  data.  These  labels  specify  the  dimensional  units  prescribed  by 
the  user  in  the  Control  Data  Block. 

If  more  than  one  page  is  required  for  an  output  group,  the  second 
and  subsequent  pages  start  with  the  page  heading  and  column  labels.  The 
title  is  not  repeated.  This  policy  is  comon  to  all  groups  of  pages  of 
printed  data. 

Table  25  represents  the  group  of  printed  output  which  lists  all  elements 
of  the  finite  element  models.  This  list  differs  from  the  printout  of  input 
blocks  by  including  elements  created  by  interior  joints.  If  an  element 
has  been  formed  because  of  a  subdivision  it  is  identified  as  "BISECTED"  and 
its  properties  are  obtained  from  the  detailed  data  given  previously  for 
the  parent  element.  If  an  element  is  not  subdivided  it  is  identified  as 
"DETAILED".  This  table  also  indicates  the  type  of  element  (reinforced 
concrete  or  wide  flange),  its  gross  length,  and  material  name. 

Table  26  Indicates  the  nature  of  the  summary  description  for  a  problem. 
This  data  defines  the  interpretation  of  input  options.  The  table  also 
indicates  additons  to  the  number  of  degrees  of  freedom  of  the  analysis 
because  of  the  characteristics  of  the  elements.  The  analyst  need  only  be 
aware  of  them  insofar  as  they  causa  reductions  in  the  storage  space  available 
for  data  of  the  mathematical  model.  In  general,  one  degree  of  freedom  is 
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CH*.»  itTFOlST  ICS  OF  TMF  SIMULtTlON 
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Table  26.  (cont'd) 
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added  for  each  frame  element.  The  frame  element  freedom  provides  for  an 
improved  deformation  model,  as  described  in  the  technical  document. 

The  self-explanatory  descriptive  summary  of  Table  26  is  the  last 
printout  of  problem  descriptive  data.  It  will  be  the  last  SINGER  printed 
output  whenever  there  are  fatal  or  critical  errors  in  the  input  data. 

If  there  are  no  errors  the  response  data  described  in  the  next  section  is 
printed. 

3.1.2  Printed  Response  Data 

These  data  describe  the  state  of  the  structure  at  particular  points 
in  time.  They  Include  groups  of  pages  which  establish  the  state  of  joint 
deflection,  which  Identify  the  internal  state  of  the  finite  elements, 
summarize  the  status  of  the  system,  and  describe  solution  progress  and 
accuracy.  The  amount  of  data  printed  in  each  group  varies  with  the  print¬ 
out  level  selected  by  the  user. 

The  minimum  number  of  load  steps  or  time  points  at  which  printouts 
are  made  is  controlled  by  the  user's  Block  1  input  data.  The  maximum 
number  depends  on  the  time  interval  used  in  numerical  integration.  In 
general,  the  state  data  is  printed  at  time  points  selected  to  attain 
the  user's  specified  relative  accuracy. 

Groups  of  printouts  occur  for  each  load  step  or  time  step.  Table 
27  lists  the  types  of  print  tables  that  may  be  provider'  in  the  order  in 
which  they  are  printed  out  in  each  printout  group.  Table  numbers  are 
those  of  this  report.  The  same  tables  are  printed  regardless  of 
whether  the  analysis  is  static  or  dynamic. 

Table  28  illustrates  printout  of  the  search  progress.  This  print¬ 
out  illustrates  the  detailed  level.  It  Includes  values  of  scaled 
variables  at  the  beginning  and  end  of  each  "linear  minimization"  -  a 
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Table  27 


TABLES  IN  A  PRINTOUT  GROUP 


Tuble 

Table 

Table 

No. 

Title 

Contenta 

28 

MINIMIZATION 

PERFORMANCE 

Valuea  of  the  work  function  and  8upportlng  data 
developed  during  aolution  search. 

30 

JOINT  STATUS 

Components  of  joint  displacements  for  all  joints. 

31 

INTERNAL 

ENERGY 

Data  defining  the  distribution  of  stored  and 
dissipated  energy  among  elements  of  the  structure. 

32,34 

AVERAGE 

STRESSES 

Data  defining  the  stress  state  in  elements  of  the 
structure. 

33,35 

SOLUTION 

ACCURACY 

Data  defining  errors  in  conservation  of  energy  and 
satisfaction  of  equilibrium. 
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Table  28 .  Minimization  Performance 
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one-dimensional  search  In  the  direction  developed  from  the  current 
measurement  of  the  gradient  and  the  estimate  of  the  curvatures  of  the 
work  function  space.  At  the  minimum  level,  this  printout  is  reduced  to 
only  values  of  the  work  function. 

Table  29  provides  more  complete  information  on  the  interpretation  of 
most  of  the  labels  used  in  the  printout  of  Table  28 .  A  special  message 
also  Indicates  the  basis  for  stopping  the  search.  This  particular  run 
indicates  that  the  run  completion  is  "NORMAL"  -  i.e.  the  minimization 
accuracy  desired  has  been  obtained.  Other  runs  may  state  "LINEAR  MINIMI¬ 
ZATION  FAILED  TO  CHANGE  THE  FUNCTION  VALUE"  or  "THE  MAXIMUM  NUMBER  OF 
LINEAR  MINIMIZATIONS  WAS  EXCEEDED."  The  former  comment  means  the  work 
function  could  not  be  changed  by  moving  in  the  last  search  direction. 

The  latter  comment  means  insufficient  minimizations  were  allowed  to 
develop  an  adequate  solution. 

Table  30  represents  the  printout  of  joint  state  data.  These  data 
define  the  state  of  motion  of  deflection  at  all  element-connected  joints 
specified  by  the  analyst.  This  printout  is  always  produced  regardless 
of  the  print  option  level.  It  is  always  the  first  response  data  printed 
after  search  completion. 

There  are  three  subgroups  of  pages  which  comprise  printed  data  of 
element  response.  These  provide  element  energy  distribution.  Internal 
stress  -strain  data,  and/or  element  Internal  forces.  Energy  is  activated 
at  the  minimum  print  level.  Stress-strain  and  internal  force  printouts 
are  activated  by  the  ISTRES  flag  of  Control  Block  input. 

Table  31  illustrates  the  element  energy  response  printout.  The 
energy  data  suggest  the  relative  Importance  of  each  element  to  the 
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Table  29 


Item 

FUNCTION  VALUE 

VARIABLES 

GRADIENTS 

CURVATURES 

FRACTION  OF 
DIRECTION  TAKEN 

FUNCTION 
SUBROUTINE  CALLS 

GRADIENT 
SUBROUTINE  CALLS 

MINIMIZATION 

ACCURACY 

TOTAL  ERROR 
MEASURES 


MINIMIZATION  PERFORMANCE  DATA 


Interpretation 


Value  (unsealed)  of  the  function  whose  minimum  is 
defined  by  the  equilibrium  configuration  of  the  system. 

Values  (translation  displacements  are  scaled  by  1.0/ 
average  element  length)  for  each  unrestrained  displace¬ 
ment  unknown. 

The  derivative  (measured  numerically)  of  the  work 
function  with  respect  to  each  of  the  variables. 

Estimates  of  the  diagonals  of  the  Influence  matrix  for 
the  system  (non-dimensional) .  These  estimates  are  exact 
after  NDF  iterations  if  the  system  is  linear. 

The  scalar  multiplier,  relative  to  that  predicted  from 
curvature  estimates,  for  each  trial  in  the  search 
direction. 

The  total  number  of  tlmea  the  total  energy  of  the  system 
has  been  evaluated  for  this  time  or  load  step  search. 


The  total  number  of  gradients  evaluated  (always  a  multiple 
of  NDF  unless  user's  step  else  guesses  were  too  small). 


The  largest  displacement  change  of  the  linear  minimization. 


The  largest  component  of  relative  energy  error  (residual 
force  times  displacement/total  energy) . 
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Table  30  .  Joint  Status  Printout 
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DX  “  displacement  in  x  direction 
DY  -  displacement  in  y  direction 
RZ  “  rotation  about  z  axis 


Table  31.  Internal  Energy  Distribution 
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behavior  of  the  system  for  the  current  time  interval.  Since  the 
solution  accuracy  is  controlled  by  controlling  the  energy  error,  these 
data  are  the  aost  meaningful  data  on  element  response.  The  recoverable 
energy  is  the  strain  energy  stored  as  elastic  work.  The  dissipated 
energy  is  that  energy  lost  due  to  yielding.  The  meaning  of  percentage 
is  clarified  by  footnotes  on  Table  31.  Data  la  printed  for  every 
element  of  the  structure.  Elements  are  listed  in  the  same  sequence  as 
In  the  element  articulation  printout  of  Table  25  . 

Tables  32  and  33  illustrate  printouts  defining  the  state  of  stress 
and  strain  and  internal  force  of  a  vide  flange  beam  element.  Tables  34  and 
35  illustrate  similar  printouts  for  a  reinforced  beam  element.  Only 
stress  resultants  are  printed  for  leaf  springs. 

Table  36  exhibits  the  solution  accuracy  summary  produced  upon 
completion  of  an  integration  step.  This  summary  is  automatically  produced 
regardless  of  the  print  level  selected.  The  data  defined  the  energy  input 
to  the  system  and  its  allocation  as  stored  and  dissipated  energy.  They 
relate  the  energy  status  and  the  change  of  status  from  that  of  the 
previously  reported  time. 

The  relative  value  of  the  energy  which  is  unaccounted  for  provides  a 
measure  of  the  cumulative  accuracy  of  the  simulation.  The  source  of  this 
error  is  inexactness  in  locating  the  true  minimum  energy  state  in  the 
solution  search  and  r<  und-of f  or  manipulation  error.  Thus,  this  error 
can  often  be  reduced  by  reducing  the  energy  tolerance  allowed  by  the  user 
in  specifying  Control  Block  input.  There  is  a  minimum  value  to  which  this 
error  measure  can  be  expected  to  be  reduced  due  to  limited  precision  in 
computer  calculations. 

In  Table  36,  residual  forces  define  the  value  of  unbalanced  force 
in  each  equilibrium  equation.  These  are  identified  with  each  joint  degree 
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’esldual  energy^total  energy 


of  freedom.  The  residual  energy  Is  the  external  work  of  the  residual 
forces  (residual  force  times  displacement  changes  in  this  step) . 

3.1.3  Messages 

The  three  types  of  printed  messages  are  described  and  illustrated  in 
Table  37.  Logic  producing  these  messages  occurs  in  practically  every 
subroutine  and  printouts  are  produced  as  this  logic  is  encountered. 

Each  error  message  indict  c  iu  parentheses  the  subroutine  name  of  its 
origin. 

Table  18  through  36  illustrate  printout  output  when  no  error  messages 
arise.  Since  error  messages  are  printed  as  they  arise,  the  pattern  of 
each  of  these  printouts  is  changed  when  errors  occur. 

3.2  RETRIEVAL  FILE  DATA 

Two  types  of  data  files  are  transferred  to  peripheral  storage  units: 
a  plot  file  and  a  calculation  recovery  file.  The  plot  file  stores  data, 
defining  the  state  of  the  structure,  for  later  printing  or  plotting. 

The  recovery  file  provides  data  for  continuing  integration  of  the  equations 
of  motion  in  another  computer  run. 

Plot  files  are  written  when  Control  Block  input  so  specified.  A 
file  is  transferred  to  the  specified  peripheral  unit  at  every  time  point 
for  which  there  is  a  printout  of  response  state  data.  This  file  contains 
a  complete  record  of  the  state  of  the  structure  at  the  end  of  each  time 
interval . 

A  recovery  file  is  written  just  before  terminating  calculations,  as 
long  as  a  recovery  peripheral  unit  has  been  specified,  and  no  fatal  or 
critical  errors  have  occurred.  This  file  is  transferred  when  the  integration 
is  completed  for  the  time  Interval  of  interest.  When  calculations  must  be 
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*Under lined  terms  vary  with  problem. 


terminated  before  integration  Is  completed,  a  recovery  file  is  also 
automatically  transferred  to  facilitate  continuation  of  calculations 
on  another  computer  run. 

The  recovery  file  contains  all  data  necessary  for  continuing 
integration  from  the  last  time  point  completed  successfully.  For 
recovery  only,  the  data  for  Input  Block  1  needs  to  be  included  in  input . 
Data  from  this  block  reinitiates  the  simulation.  Accordingly,  input 
for  this  block  must  be  changed  to  reflect  the  continuation  conditions 
as  described  by  the  footnotes  of  Table  2. 

Table  38  provides  a  summary  of  all  input  data  blocks  and  groups 
of  output  data.  This  table  identifies  the  data  saved  in  the  plot  and 
recovery  files.  Retrieval  files  contain  only  binary  coded  data.  Each 
of  these  files  is  written  by  a  single  FORTRAN  WRITE  statement. 


Table  38 


RETRIEVAL  FILE  DATA 


Data 

Plot  File 

Recovery  File 

Problem  Descriptive  Data 

Control  Block  Data 

No 

No 

Remaining  Input  Blocks 

No 

Yes 

Problem  summary  data 

No 

Yes 

Joint  Deflection  Data 

Deflection,  and  velocity 
in  each  degree  of  freedom 

Yes 

Yes 

Acceleration  in  each  degree 
of  freedom 

Yes 

No 

Element  State  Data 

Energy  distribution 

Yes 

No 

Internal  loads 

Yes 

No 

Stresses  and  strains 

Yes 

Yes 

Element  status  summary 

No 

No 

Solutions  Progress  Data 

Summary  data 

No 

Yes 

Joint  balance  data 

No 

No 

Solution  progress  data 

No 

No 

Error  Messages 

No 

No 

SECTION  4 


PROBLEM  SETUP 

The  computer  running  tine  and  analysis  efficiency  depend  on  how  the 
user  describes  the  problem  using  the  input  data.  The  purpose  of  this 
section  is  to  give  a  user  some  general  guidelines  for  setting  up  the 
problem  and  controlling  the  solution  accuracy. 

In  preparing  input,  one  of  the  user's  objectives  is  to  define  the 
problem  for  the  computer  so  that  it  can  develop  accurate  results  effi¬ 
ciently.  Efficiency  is  especially  important  in  nonlinear  analysis  where 
excessively  detailed  models  can  result  in  an  increase  of  orders  of  magni¬ 
tude  in  processing  time.  These  excesses  can  be  reduced  by  using  the 
analysts'  understanding  of  the  program  and  the  structural  behavior  to  simplify 
the  problem  definition. 

In  decreasing  order  of  importance  to  processing  time,  the  user's 
objective  should  be  to  minimize  the  number  of  degrees  of  freedom,  the 
complexity  of  element  representation,  and  the  solution  search  time. 

4.1  REDUCING  THE  NUMBER  OF  EQUATIONS 

For  all  but  the  smallest  problems,  the  majority  of  the  computer 
resources  are  devoted  to  solving  the  nonlinear  equations.  Since  the  number 
of  calculations  for  solution  is  a  function  of  the  number  of  equations 
(experience  suggests  a  factor  of  two  to  ten) ,  reducing  the  number  of 
equations  is  important. 

Reducing  equations  means  modeling  with  the  fewest  joints  and  degrees 
of  freedom  per  joint.  Where  app' icr>ble,  only  half,  or  a  quadrant  of  the 
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system  should  be  modeled  using  symmetry  or  assycmetry  boundary  conditions 
to  imply  the  rest  of  the  system.  Joints  which  will  not  move  significantly 
in  a  given  direction,  such  as  the  axial  motions  of  a  laterally-loaded 
continuous  beam,  should  be  restrained.  Distributed  loadings  should  be 
chosen  rather  than  introducing  interior  joints  for  a  centrally-loaded 
span,  if  the  central  load  work  will  be  comparable  to  the  work  of  statically 
equivalent  end  loads.  Joints  that  are  close  together,  compared  with  other 
joints,  should  be  combined  in  a  single  joint.  If  the  structure  becomes 
broken  into  two  parts,  each  part  should  be  separated  to  reduce  total 
calculations. 

4.2  •  REDUCING  THE  COMPLEXITY  OF  ELEMENT  REPRESENTATION 

The  number  of  calculations  in  equation  solving  is  directly  propor¬ 
tional  to  the  number  of  energy  (and  gradient)  evaluations.  Since  each  of 
these  involves  determining  the  energies  in  each  element,  reducing  the  number 
and  complexity  of  element  models  directly  reduces  data  processing  time. 

Several  modeling  devices  are  available  with  SINGER  for  reducing  element 
complexity.  Since,  in  many  impact  problems,  much  of  the  system  behaves  in 
the  linear  stress-strain  range  throughout  the  loading  history,  the  user  is 
permitted  to  designate  that  energy  calculations  can  be  simplified  to  the 
linear  range  for  particular  elements.  Judicious  use  of  leaf  springs,  in¬ 
stead  of  modeling  with  other  line  elements,  replaces  models  which  require 
Integration  over  the  volume  of  the  element,  with  a  matrix  multiplication 
involving  relatively  few  calculations.  Similarly,  grouping  longitudinal 
steel  into  fewer  groups  and  using  videflange  beams  Instead  of  reinforced 
concrete  models  reduces  energy  Integration  calculations. 
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4.3  REDUCING  THE  SOLUTION  SEARCH  TIME 

The  solution  search  tine  is  measured  by  the  number  of  linear  minimi¬ 
zations  needed  to  arrive  at  an  acceptable  solution  of  the  equations  of 
equilibrium  guidelines  for  modeling  to  minimize  the  time  center  around 
providing  good  guesses  for  the  load  or  time  step  and  reducing  the  complexi¬ 
ty  of  nonlinearity. 

4.3.1  Selection  of  Load  or  Time  Steps 

In  obtaining  response  to  static  loadings,  the  load  can  be  applied  over 
a  number  of  steps,  If  physically  meaningful.  This  approach  tends  to  reduce 
search  by  keeping  configuration  changes  small  enough  so  convergence  Is 
rapid  for  each  step  and  thereby,  the  total  search  time  is  reduced.  If 
load  steps  are  too  large,  the  user  will  observe  that  progress  In  reducing 
the  work  function  may  be  slow.  This  difficulty  is  often  eliminated  simply 
by  restarting  the  analysis  search  using  the  recovery  feature  of  SINGER. 

For  dynamic  loadings,  search  time  Is  reduced  by  requiring  Integration 
time  steps  of  the  order  of  one-fifth  the  shortest  period  of  resonance. 

This  eliminates  needless  time  cutting  and  multiplying  that  will  occur  If 
SINGER  determines  the  appropriate  time  step  by  trial  and  error.  Some  of 
this  trial  and  error  will  occur  naturally  when  the  structure  goes  plastic. 
Then,  the  search  can  be  reduced  by  restarting  with  user  supplied  time 
steps. 

4.3.2  Reducing  the  Complexity  of  Nonlinearity 

The  complexity  of  nonlinearity  Is  reduced  by  limiting  the  type  of 
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nonlinearity,  simplifying  material  representations,  and,  in  dynamic  response 
analysis,  by  simplifying  the  mass  model. 

Types  of  nonlinearity  include  orientation  and  material  nonlinearities. 
Designating  that  only  small  angle  changes  can  occur  at  joints  (sine  (angle)  » 
angle,  cosine  (angle)  *  1.0)  reduces  the  equations  to  representing  nonlinear 
strain-displacement  terms  in  geometry.  Since  in  many  problems  these  terms 
only  weakly  couple  axial  and  lateral  motions,  the  equations  become  nearly 
linear  and  the  solution  search  is  achieved  in  less  than  NDF  search  steps 
for  each  new  loading,  where  NDF  is  the  number  of  degrees  of  freedom. 

Material  nonlinearities  are  eliminated  by  modeling  each  material  with 
a  single  line.  Since  folded  line  models  involve  a  jump  discontinuity  in 
effective  Young's  Modulus  at  each  fold,  the  single  line  eliminates  approxi¬ 
mations  made  over  the  discontinuity  and  consequent  added  search  steps. 
Likewise,  extending  the  domain  of  material  linearity,  reducing  the  number 
of  folds,  and  minimizing  the  fold  angle  change  reduce  the  number  of  search 
steps. 

Simplifying  the  mass  model  usually  involves  reducing  the  total  wither 
of  masses  by  lumping  at  fewer  joints.  The  effect  of  this  reduction  is  to 
Increase  the  shortest  resonance  period  and  consequently  permit  larger  time 
steps  in  dynamic  response  prediction.  Since  the  shortest  period  can  be 
bounded  by  the  period  associated  with  movement  of  only  the  smallest  mass, 
reducing  the  number  of  masses  should  usually  involve  increasing  the  mag¬ 
nitude  of  thi  smallest  mass. 


88 


